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ABSTRACT:  The  friction  coefficient  and  velocity  profile  are 

calculated  fron  the  relation  between  the  shear,  the  sun  of  the 
Molecular  and  eddy  viscosity,  and  the  velocity  gradient.  The 
non-dlmenslonal  shear  distribution  Is  assuaed  to  be  fixed.  The 
eddy  viscosity  across  the  boundary  layer  Is  obtained  by  Joining 
a  distribution  for  the  wall  region  to  one  for  the  outer  region. 
By  use  of  the  shear  and  eddy  viscosity  distributions  the  veloc¬ 
ity  profile  Is  calculated  from  the  wall  to  the  outer  edge  of  the 
boundary  layer  for  all  Reynolds  numbers  without  using  the  con¬ 
cept  of  laminar  sublayer,  transition  region,  logarithmic  region, 
etc . 

The  heat  and  mass  transfer  coefficients  and  the  temperature  and 
concentration  profiles  are  calculated  by  a  similar  method. 
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SUMMARY 

A  new  method  for  the  calculation  of  mean  flow  quantities  In 
a  constant  property  turbulent  boundary  layer  on  a  flat  plate  is 
developed.  In  this  method  the  friction  coefficient  and  velocity 
profile  are  calculated  from  the  relation  between  the  shear,  the 
sum  of  the  molecular  and  eddy  viscosity,  and  the  velocity  gradient. 
The  shear  distribution,  experimentally  obtained  by  previous 
investigators,  is  assumed  to  be  fixed.  The  eddy  viscosity  is 
obtained  by  Joining  a  distribution  for  the  wall  region  to  one 
for  the  outer  region  of  the  boundary  layer.  Both  distributions 
were  calculated  by  previous  investigators  from  measured  flow 
quantities.  By  use  of  the  shear  and  eddy  viscosity  distributions 
the  friction  coefficient  is  calculated  and  agrees  well  with 
accepted  values  over  the  entire  range  of  Reynolds  number.  For 
large  Reynolds  numbers  the  logarithmic  friction  formula  is 
obtained. 

The  velocity  profile  is  calculated  by  the  same  method  as 
the  friction  coefficient.  The  entire  profile,  from  the  wall  to 
the  outer  edge,  is  calculated  without  using  the  concept  of 
laminar  sublayer,  transition  region,  logarithmic  region,, etc. 

The  calculated  velocity  profile  changes  gradually  from  the 
laminar  type  at  small  Reynolds  numbers  to  the  turbulent  type  at 
larger  Reynolds  nxuabers.  For  sufficiently  large  Reynolds  numbers 
the  logarithmic  formula  is  obtained  for  one  part  of  the  velocity 
profile  and  the  velocity  defect  formula  for  another  part. 

The  heat-transfer  coefficient  and  non-dimensional  tempera¬ 
ture  profile  for  small  temperature  differences  are  calculated 
by  the  same  method  as  the  friction  coefficient  and  velocity 
profile.  In  these  calculations  the  non-dimensional  heat-transfer 
distribution,  which  roust  be  known,  is  assiimed  to  be  the  same  as  the 
non-dimensional  shear  distribution.  In  addition  to  the  Prandtl 
number,  the  ratio  of  the  eddy  thermal  dlffusivlty  to  the  eddy 
kinematic  viscosity  also  appears.  This  ratio  is  calculated  from 
experimental  data  obtained  by  previous  investigators.  The  calcu¬ 
lated  ratio  of  local  Stanton  number  to  half  the  local  friction 
coefficient  is  found  to  vary  from  about  1.22  at  a  Reynolds  number 
of  about  3.5  X  10^,  based  on  the  distance  to  the  leading  edge, 
to  about  1.09  at  a  Reynolds  number  of  2.2  x  10^.  The  calculated 
non-dimensional  temperature  profiles  differ  slightly  from  the 
velocity  profiles. 

Turbulent  flow  over  a  surface  on  which  condensation  or 
evaporation  is  occurring  is  treated  by  assuming  that  the  rate  of 
condensation  or  evaporation  is  small,  that  the  Schmidt  number  is 
equal  to  the  Prandtl  number,  and  that  the  eddy  dlffusivlty  is 
equal  to  the  eddy  thermal  dlffusivlty.  The  mass-transfer  coef¬ 
ficient  is  then  equal  to  the  heat- transfer  coefficient,  and  the 
non-dimensional  concentration  profile  is  the  same  as  the  non- 
dimensional  temperature  profile. 
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INTRODUCTION 

The  problem  of  calculating  the  friction  coefficient  and 
velocity  distribution  for  a  constant  property  turbulent  boundary 
layer  is  an  old  one.  In  a  rigorous  sense  the  problem  is 
unsolved  and  will  remain  so  until  shear  turbulence  is  understood 
better.  At  present  the  friction  coefficient  is  usually  calculated 
from  a  formula  designed  to  fit  either  experimental  data  or  the 
logarithmic  formula  (refs.  (1)  and  (2)).  The  logarithmic 
formula  is  theoretically  derived  but  contains  two  constants  that 
must  be  found  from  experiment.  This  formula  can  be  derived  in  a 
number  of  ways,  for  extunple,  by  use  of  Prandtl's  or  von  Karman's 
relations  for  the  eddy  viscosity  or  by  Millikan's  method  of 
overlap  (ref.  (1)). 

The  analyses  that  lead  to  a  logarithmic  friction  formula 
also  result  in  the  logarithmic  velocity  profile.  This  velocity 
distribution  agrees  reasonably  well  with  experiment  over  a  portion 
of  the  boundary  layer  from  near  the  edge  of  the  laminar  sublayer 
to  about  one-fifth  of  the  boundary  layer  thickness  from  the  wall. 
The  region  very  close  to  the  wall  requires  a  separate  treatment 
that  gives  the  result  that  In  this  region  the  velocity  is 
directly  proportional  to  the  distance  from  the  wall.  The 
velocity  profile  must  vary  continuously  from  the  linear  to  the 
logarithmic,  but  the  variation  is  not  obtainable  from  either  the 
procedure  that  gives  the  linear  profile  or  from  that  which 
yields  the  logarithmic  profile.  A  number  of  Investigators  have 
treated  this  problem  (ref.  (1)).  To  calculate  the  velocity 
distribution  thus  requires  that  the  boundary  layer  be  divided 
into  a  number  of  regions.  This  is  true  even  if  Cole’s  correla¬ 
tion  of  velocity  profiles  (ref.  (1))  Is  used.  A  complete 
review  and  detailed  discussion  of  the  present  knowledge  of  the 
constant  property  turbulent  boundary  layer  on  a  flat  plate  is 
given  in  reference  (1). 

In  the  present  analysis  the  velocity  profile  and  friction 
coefficient  are  obtained  from  the  relation  between  the  shear, 
the  sura  of  the  molecular  and  eddy  viscosity,  and  the  velocity 
gradient.  The  shear  distribution,  experimentally  obtained  by 
previous  investigators,  is  assumed  to  be  fixed.  The  eddy 
viscosity  is  obtained  by  Joining  a  distribution  for  the  wall 
region  to  one  for  the  outer  region  of  the  boundary  layer.  Both 
distributions  were  calculated  by  previous  investigators  from 
measured  flow  quantities.  By  use  of  the  shear  and  eddy  viscosity 
distributions  the  friction  coefficient  is  calculated  and  agrees 
well  with  accepted  values  over  the  entire  range  of  Reynolds 
number.  For  the  limiting  case  of  imro  Reynolds  number,  the 
friction  coefficient  becomes  almost  equal  to  the  exact  Blasius 
value  for  laminar  flow  (ref.  (2)).  At  the  other  extreme,  very 
large  Reynolds  number,  the  analytic  expression  for  the  logarithmic 
friction  formula  is  obtained. 
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The  velocity  profile  is  calculated  by  the  same  method  as 
the  friction  coefficient.  Unlike  the  usual  procedure,  the  entire 
profile  from  the  wall  to  the  outer  edge  Is  calculated  without 
splitting  the  boundary  layer  in  a  more  or  less  arbitrary  manner 
into  laminar  sublayer,  transition  region,  logarithmic  region, 
etc.  The  calculated  velocity  profile  changes  gradually  from 
almost  the  exact  Blaslus  laminar  profile  in  the  limit  of  zero 
Reynolds  number  to  the  turbulent  type  at  larger  Reynolds  numbers. 

The  heat-transf er  coefficient  and  temperature  profile  lor 
small  temperature  differences  are  calculated  in  a  manner  similar 
to  that  for  the  friction  coefficient  and  velocity  profile.  That 
is,  Instead  of  the  relation  between  the  local  shear,  the  total 
viscosity,  and  the  velocity  gradient,  there  Is  used  the  relation 
between  the  local  heat  transfer,  the  sum  of  the  molecular  and 
eddy  conductivity,  and  the  temperature  gradient.  The  information 
required  to  calculate  the  eddy  conductivity  is  obtained  from  the 
experimental  data  of  reference  (3) .  The  non-dimensional  heat- 
transfer  distribution  across  the  boundary  layer  is  assumed  to 
be  the  same  as  the  non-dimensional  shear  distribution. 

For  small  concentration  of  a  foreign  gas,  for  Schmidt 
number  equal  to  Prandtl  number,  for  eddy  dlffuslvity  equal  to 
the  eddy  thermal  dlffuslvity,  and  for  the  non-dimensional  mass- 
transfer  distribution  across  the  boundary  layer  equal  to  the 
non-dimensional  heat- transfer  distribution,  the  concentration 
profile  and  mass-transfer  coefficient  can  be  obtained  directly 
from  the  temperature  profile  and  heat-transfer  coefficient. 


ANALYSIS 

Friction  Coefficient  and  Velocity  Profile 

In  laminar  boundary  layer  flow  the  shear  stress  is  related 
to  the  velocity  gradient  by  the  relation 

•7*  —  lA  ^ 

For  turbulent  flow,  it  is  assumed  that  an  eddy  kinematic 
viscosity  Cg)  can  be  found  such  that 

(2) 

Equation  (2)  can  be  written  as 

"  -C  5-  ■ 
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or, 


for 


constant  property 
T  _ 


^0. 


(3) 


From  (3) 


it 


follows  that 


o 


Therefore 


(4) 


(5) 


T 

To  evaluate  the  integral  in  (5) , 


and  must  be  known 


from  the  wall  to  the  outer  edge  of  the  boundary  layer.  The 
ratio  ,  shown  In  figure  (1) ,  was  obtained  from  figure  7-13 

of  reference  (1);  figure  7-13  gives  the  ratio  of  the  measured 
turbulence  sh^r  stress  to  the  wall  shear  stress.  Near  n  -  0, 
the  curve  of  was  drawn  slightly  above  the  curve  of  figure 

W 

7-13  of  reference  (1)  in  order  to  agree  more  closely  with  the 
experimental  points  shown  in  that  figure  and  also  to  allow  for 
the  viscous  shear  and  to  give  the  slope  its  required  value  of 
zero  at  the  wall.  In  figure  1  is  also  shown  the  ratio  ^  cal- 

culated  by  Fediaevsky's  method  (ref.  (4)).  In  this  method  the 
ratio  ^  is  expressed  as  a  fourth  degree  polynomial  in  ti  and 

the  coefficients,  which  depend  on  the  value  of  i  and  its 

derivatives  at  n  —  0  and  at  n  -  1,  are  evaluated  by  use  of  the 
boundary  layer  equation  of  motion.  Although  Fediaevsky's 
distribution  of  ^  is  not  very  far  from  the  measured  one,  the 

measured  one  was  used  in  the  present  analysis.  For  a  flat  plate, 
Fediaevsky's  method  makes  a  function  of  n  alone;  the  same 

assumption  is  made  in  this  analysis. 


The  ratio  IT  is  obtained  from  reference  (1).  The  discus¬ 
sion  on  page  526  and  the  data  of  figure  7-41  of  reference  (1) 
indicate  that  very  near  the  wall 
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Relation  (6)  is  universal  under  the  assumption  that  the  portion 
of  the  velocity  boundary  layer  very  near  the  wall  can  be  com¬ 
pletely  described  by  the  quantities  u,  em,  p,  \i ,  r^,  and  y. 

Further  from  the  wall,  the  discussion  on  page  493  and  the 
data  In  figure  7-17  of  reference  (1)  Indicate  that 


-b 


Because  this  relation  does  not  contain  the  viscosity.  It 
probably  applies  only  at  large  Reynolds  numbers.  At  small 
Reynolds  numbers,  Is  probably  a  function  of  Re^  In  addition 

to  ri,  u.r 

Because  (6)  applies  only  near  the  wall  and  (7)  applies 
only  further  out,  (5)  Is  written  as 


T_ 

J 


where  Is  the  smallest  value  of  r|  for  which  (7)  applies.  For 
^  depends  on  y.  Instead  of  on  r|.  Consequently,  the 

first  Integral  In  (8)  Is  written  as 


where 


7?e  -  ^ 


j*-  ■“ 


_  _ 


=  1 


'■J*,!  = 
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The  notation  ~  points  out  that  is  assumed  to  be  a 

function  of  near  the  wall  as  well  as  further  out. 

Equation  (8)  then  becomes 

°  1, 

The  function  f2(ri)  was  obtained  directly  frcMn  the  faired 
curve  in  flf^ure  7-17  of  reference  (1)  and  is  shown  in  figure  2. 

The  function  fi(y+)  was  obtained  from  figures  7-16  and  7-41  of 
reference  (1)  and  is  shown  in  figure  3.  The  points  are  the  data 
points  of  flgrures  7-16  and  7-41  of  reference  (1).  For  15  <  y^.  <  52, 
the  curve  of  fj  in  figure  3  is  a  slightly  adjusted  version  of  the 
curve  of  figure  7-16  of  reference  (1) .  The  data  of  reference  (1) 
do  not  Indicate  an  upper  limit  to  the  rising  portion  of  In 

the  present  analysis,  however,  fj^  is  postulated  to  have  the 
constant  value  .393  for  a  range  of  y+  >  52  if  Ref  is  sufficiently 
large.  The  function  fj^  Is  postulated  to  have  this  constant 
portion  because  if,  at  a  sufficiently  large  value  of  Ref,  f2  is 
plotted  against  y^  by  use  of  (9)  and  the  identity 


the  portion  of  f2  ^or  n  less  than  about  ,06  in  figure  2  with  a 
slope  of  .393  lies  on  the  line  “  .393  in  figure  3.  Therefore, 

€m/u,y  would  be  double  valued  for  y^.  >  52  unless  fi  either  did 
not  extend  beyond  y^.  -  52  or  were  equal  to  f2  beyond  y+  ~  52.  By 
taking  f^  -  .393  for  y^  >  52,  it  is  not  necessary  to  limit  the 
range  of  fi  by  the  fixed  number  52.  Moreover,  making  f^  extend 
beyond  y+  -  52  necessitates  that  f2  “  Kri  in  the  region  of  the 
boundary  layer  in  which  f^  and  f2  apply  simultaneously.  Thus, 
from  (6)  it  follows  that 

Therefore,  where  (6)  and  (7)  are  both  valid  it  is  necessary  that 
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At  high  Reynolds  numbers  the  upper  limit  of  is  thus  postulated 
to  be  the  largest  value  of  y+  for  which  f^  =  K,  that  is,  at 
about  .06  Ref. 

€m 

The  lower  limit  rif  for  ±2  is  found  by  first  calculating 
from  f2(Tl)  by  (11)  and  then,  by  use  of  (9),  plotting  these  * 
values  of  together  with  ff  as  shown  in  figure  4.  The 

extension  of  f2  down  to  y+  .=  0  is  merely  for  the  purpose  of 
calculating  .  When  f^  and  f2  intersect  at  y^  <  52,  the  value 
of  y+  at  the  intersection  is  y^.  f.  Equation  (9)  then  gives  ni . 
When  there  is  no  intersection  at  y+  <  52,  y+,l  is  52  and  is 

.  The  curve  of  f2  for  Ref  =  2960  in  figure  4  is  an  example 

of  a  Reynolds  number  that  is  so  large  that  ff  and  f2  do  not 
intersect  for  y+  <  52.  Therefore,  52  is  used  for  y+,1.  The 
curve  for  Ref  -  200  shows  the  behavior  at  a  low  Reynolds  number. 
Here,  ff  and  f2  intersect  at  y+  <  52.  When  the  Reynolds  numbers 
are  so  low  that  y+,i  <  52,  then  y+,i  depends  on  Ref. 


By  use  of  figures  1,  2,  and  3,  and  this  method  for  finding 
y+  1  and  the  value  of  Cf/2  was  calculated  from  (10)  for  six 

values  of  Ref.  The  values  of  y^.  i  and  and  the  corresponding 
values  of  fi  and  f2  for  the  three  values  of  Ref  for  which  y+  i 
is  less  than  52  are  given  in  Table  I.  In  Table  II  are  listed 
the  values  of  Cf/2  and  Re5.  The  value  of  Cf/z  for  Ref  -  2960, 
namely  .037S,  is  almost  the  same  as  the  value  .037  given  in 
reference  (1).  In  figure  5  the  calculated  values  of  Cf/2  are 
compared  with  those  obtained  by  use  of  the  logarithmic  velocity 
profile  and  given  in  Table  21.1  of  reference  (2).  The  agreement 
is  good  over  the  entire  range  of  Ref.  It  is  noted,  however, 
that  because  (7)  does  not  contain  the  viscosity,  it  is  probably 
valid  only  at  large  Reynolds  numbers.  At  small  Reynolds  numbers 
the  function  f2  Is  probably  somewhat  different  from  that  shown 
in  figure  2  with  the  result  that  Cf/2  is  also  probably  different 
from  that  calculated. 


The  velocity  profile  is  obtained  by  solving  equation  (3) 
for  and  integrating  between  0  and  r.  Thus, 


cl»|^  (13) 
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For  the  same  reason  that  (5)  was  written  as  (8),  (13)  is  written 


for  y^,  <  y+,1,  and  as  ^ 

'1'  “  /^T+Y^  -f-  ■%  /  I  1  (15) 

for  T|  >  T|i .  The  velocity  profiles,  4^(r))  ,  are  tabulated  in  Table 
III  and  shown  in  figure  6  for  the  six  values  of  Ref  for  which 
Cf/2  was  calculated.  Because  there  is  a  unique  relation  between 
Ref  and  Cf/2,  it  follows  that  the  velocity  profile  is  uniquely 
related  to  Cf/2.  Included  in  Table  III  and  figure  6  is  the 
profile  for  laminar  flow.  This  profile  was  calculated  by  using 
the  same  IL-  distribution  as  for  turbulent  flow,  that  shown  in 

figure  1,  and  by  putting  equal  to  zero  in  equation  (5)  to  get 

This  result  was  then  used  in  (13)  to  get  _ 

ru  r\ 


626 


The  velocity  profiles  in  figure  6  show  a  gradual  change  from  the 
laminar  to  the  turbulent  type  as  Ref  Increases. 

Note  that  both  the  friction  coeffi-ient  and  velocity  profile 
become  exactly  those  for  laminar  flow  when  the  Reynolds  number 
becomes  zero.  Thus,  fa  in  (5)  and  (13)  can  be  written  as_ 

/  )  ,|“cx  *  1?  ^  •  Then,  as  Re*  decreases  to  zero  £:?  also 

(uvS'JlI'^r  r-^ 

decreases  to  zero  because  the  friction  Reynolds  number,  'llCf/2  Re5, 
is  zero  for  Re^  equal  to  zero.  Equations  (5)  and  (13)  then 
become  (16)  and  (17),  respectively,  the  equations  for  laminar 
flow . 

In  figure  7  all  the  velocity  profiles  of  figure  6  except 
the  laminar  one  are  shown  in  the  form  u+  against  y^.  The  relation 
between  and  ^  is 

4 
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and  the  relation  between  and  n  is  given  by  (9) .  Each  velocity 
profile  has  a  straight  line  portion  that  lies  on  the  single  line 
which  is  fitted  by  the  equation 

U.4_=  Z.-5A  \n  y_^  -V-  5.8  9 

Each  velocity  profile  leaves  this  line  at  a  sufficiently  large 
value  of  u+.  Calculations  indicate  that  this  occurs  when  cb 
exceeds  about  .75.  ' 


The  profile  for  Re^  =  25  lies  markedly  below  all  the  others. 
The  reason  is  that  25  is  so  small  a  value  of  Ref  that  large 
values  of  r  occur  for  small  values  of  y+  (see  eq.  (9)).  Conse¬ 
quently,  ^  -  ,  which  is  assumed  to  depend  on  alone,  is  smaller 

over  the  range  of  integration  in  (15)  than  it  is  for  larger  values 
of  Ref.  For  example,  at  y+  =*  12.2,  which  is  y+,  i  in  this  case, 
n  is  large  enough  for  to  equal  .600,  a  value  appreciably  less 

than  unity.  As  a  result  U4.  is  small  over  the  entire  boundary 
layer  and  its  maximum  value  is  only  9.31. 

In  figure  7  are  also  shown  Klebanoff's  data  for  Ref  -  2960 
from  figure  21  of  reference  (5).  Klebanoff’s  data  differ 
slightly  from  the  curve  of  the  present  analysis  for  values  of 
y+  between  about  20  and  120.  If  the  curve  of  ff  of  reference 
(1)  had  been  used  without  refairing  it,  better  agreement  would 
have  been  obtained  with  Klebanoff's  data. 


In  figure  i  are  shown  the  computed  velocity  profiles  in  the 

form  We— against  n.  Also  shown  is  the  curve  from  figure  7-4  of 

■ 

reference  (1);  this  curve  is  a  mean  through  experimental  data. 
There  is  fair  agreement  between  the  calculated  velocity  profiles 
and  this  mean  curve  for  the  three  largest  values  of  Ref  but  not 
for  the  three  smallest  values.  The  three  smallest  Reynolds 
numbers  are  probably  much  smaller  than  those  for  which  data  were 
used  to  construct  the  curve  of  figure  7-4  of  reference  (1).  The 
comparisons  with  experiment  in  figures  7  and  i  indicate  that  a 
calculated  profile  probably  differs  slightly  from  an  experimental 
one  at  the  same  friction  Reynolds  number. 


and 


The 


From 

were 


ratio 


the  computed  velocity  profiles  the  quantities 
calculated  numerically  from  their  definitions 


was  calculated  from  — 

e 


B 

,  namely, 
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These  quantities  are  given  in  Table  IV.  By  use  of  ,  the 

value  of  Re^,  a  more  convenient  quantity  than  Re^  or  Ref,  was 
calculated  for  each  value  of  Ref  for  which  Cf/2  was  computed. 
These  values  of  Req  are  given  in  Table  IV  and  the  dependence  of 
Cf/2  on  Reo  is  shown  in  figure  9. 


When  Cf/2  is  known  as  a  function  of  Req,  the  dependence  of 
Cf/2  and  of  Cp  on  Re^  can  be  calculated.  The  procedure  is  to 
use  the  boundary  layer  momentum  equation  in  the  form 


d  “f?#  X  ^ 


(19) 


which,  upon  integration,  gives 


Ce=0  at  x=o) 


(20) 


The  relation 


is  also  used. 


'Re# 


(21) 


To  calculate  Rex  from  (20),  Cf/2  must  be  known  as  a  function 
of  Rea  for  small  values  of  Rea*  In  the  present  analysis  the 
turbulent  boundary  layer  gradually  becomes  laminar  as  Rsq 
approaches  zero  so  that  it  is  unnecessary  to  assume  either  a 
turbulent  boundary  layer  at  the  leading  edge  or  a  transition 
position.  In  order  to  use  (20),  Cf/2  was  extrapolated  from  its 
value  at  Req  equal  to  27.7  down  to  zero  by  fitting  the  equation 


to  the  values  of  Cf/2  at  Req  >  27.7  and  191  (see  Tables  II  and 
IV).  It  was  found  that  ko  ••  .0789  and  n^  ■■  .5785.  For  Req  ••  100 
equation  (20)  then  gives  Rex  ~  11,530  and  equation  (21)  gives 
Cp  -  .01735.  Although  kg  and  ng  should  aporoach  their  laminar 
values,  namely,  .2205  and  1,  instead  of  remaining  fixed  as  Req 
becomes  zero,  it  can  be  shown  that  at  Req  ~  100  the  error  in  Rex 
caused  by  keeping  kg  and  ng  fixed  is  no  larger  than  2  percent. 

The  percentage  error  in  Rex  and  Cp  decreases  as  Req  Increases. 

For  values  of  Req  larger  than  100,  equation  (20)  was  Integrated 
numerically  by  use  of  figure  9.  In  figure  10  is  shown  the 
calculated  dependence  of  Cp  on  Rex;  also  shown  is  the  dependence 
of  Cp  on  Rex  predicted  by  the  Prandtl-Schllchtlng  and  Schultz- 
Grunow  formulas  (ref.  (2)).  The  agreement  is  good. 
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la  figure  11  is  shown  the  connection  between  Roq  and  Rex 
calculated  by  use  of  (20).  Figure  11  is  given  because  it 
enables  quantities  known  as  a  function  of  Req  to  be  easily 
obtained  as  a  function  of  Rex  and  vice  versa.  Thus,  from  this 
figure  and  figure  9  the  dependence  of  Cf/2  on  Rex  can  be  found. 


Heat-Transfer  Coefficient  and  Temperature  Profile 


The  heat-transfer  coefficient  and  temperature  profile  are 
calculated  in  the  sane  way  as  the  friction  coefficient  and 
velocity  profile.  The  expression  for  the  heat  transfer  by 
conduction  in  a  direction  normal  to  the  wall  in  a  laminar  flow 
is 


H  = 


(22) 


It  Is  assumed  that  an  eddy  thermal  diffusivity  exists  that  allows 
the  heat  transfer  in  a  turbulent  flow  to  be  written  as 

i  =  (23) 

Equation  (23)  can  be  written  as 


After  defining  a  non-dimensional  temperature  t  by 

Te  -Tw  > 

introducing  a  temperature  boundary- layer  thickness  5'p,  and 
dividing  by  q^,  (24)  becomes 


q  _  K  (Te^'v) 


vV  L-t 


(25) 


T 


where  nr.  -  .  After  Introduction  of  the  Stanton  number 

C  =  _  V  ^  , 

equation  (25)  becomes  ' 
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From  (26)  it  follows  that 


=  / 


Before  (28)  can  be  used  to  calculate  the  ratio  must 

be  known  as  a  function  of  np.  Unlike  -2-  ,  no  measurements  of 

q  i  iw  ’ 

-3-  on  a  plate  could  be  found.  Consequently,  to  obtain  an 
approximate  distribution  of  ,  the  Fediaevsky  method  for 

getting  was  used  to  investigate  the  dependence  of  on  np 

This  method  gave  a  reasonable  approximation  to  ^  (see  fig.  1) 

Two  conditions  to  evaluate  the  fourth  degree  polynomial  of 
Fediaevsky 's  method  are 

^=1  y  =  o 

and  . 


-fc  = 


Also,  from  the  energy  equation, 


=-21- 


and  the  conditions,  Vy„  -  0  and 


'ft)  - 


(II-).-  o. 


it  follows  that  both 


Moreover,  f or/^^  \  —  0,  the  condition  at  y  —  6't>  is 

^  ft)  =-  .  ^ 

These  five  conditions  on  exactly  the  same  as  on  Tw 

Consequently,  if  2:  and  _1.  are  approximated  by  fourth  degree 

q  Tw  "w  -p 

polynomials  is  exactly  the  same  function  of  *7p  as  is  of 

n.  Therefore,  the  relation  between  and  is  assumed  to  he 

the  same  as  the  experimental  one  between  "tQ  and  p.  Actually, 
however,  -3-  depends  slightly  on  the  Prandtl  number  (see  fig. 
7-55,  rof.^'d)). 
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In  figure  11  is  shown  the  connection  between  Rsq  and  Rex 
calculated  by  use  of  (20).  Figure  11  is  given  because  it 
enables  quantities  known  as  a  function  of  Req  to  be  easily 
obtained  as  a  function  of  Rex  and  vice  versa.  Thus,  from  this 
figure  and  figure  9  the  dependence  of  Cf/2  on  Rex  can  be  found. 

Heat-Transfer  Coefficient  and  Temperature  Profile 

The  heat-transfer  coefficient  and  temperature  profile  are 
calculated  in  the  same  way  as  the  friction  coefficient  and 
velocity  profile.  The  expression  for  the  heat  transfer  by 
conduction  in  a  direction  normal  to  the  wall  in  a  laminar  flow 
is 

S  =  -  K  1“  (22) 

It  is  assumed  that  an  eddy  thermal  dlffusivity  exists  that  allows 
the  heat  transfer  in  a  turbulent  flow  to  be  written  as 

(23) 

Equation  (23)  can  be  written  as 


After  defining  a  non-dimensional  temperature  t  by 

.  ^  ~Tw 

Te  -Tvv  > 

introducing  a  temperature  boundary- layer  thickness  6'f,  and 
dividing  by  q^,  (24)  becomes 


wh.re  .  After  Introduction  of  the  Stanton  number 

S  =  -  V  , 

*  P^d'T.CTe-Te/') 

equation  (25)  becomes  ' 


_3_  _  _J _ ; —  f  I  4-  — 

p  U  Vv;UJ  J  -i'll 
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From  (26) 


it  follows  that 

'_3- 

/  sv\-  1 


S's.  S' 


ci*lT 


(27) 


or 


A 


(28) 


Before  (28)  can  be  used  to  calculate  the  ratio 


■fe 


must 

be  known  as  a  function  of  nr.  Unlike  ,  no  measurements  of 

Q  ‘w  ’ 

i  on  a  plate  could  be  found.  Consequently,  to  obtain  an 
iw  q 

approximate  distribution  of  q  ,  the  Fediaevsky  method  for 
X  a 

getting  was  used  to  investigate  the  dependence  of  on  np 

This  method  gave  a  reasonable  approximation  to  ^  (see  fig.  1) 

‘W 

Two  conditions  to  evaluate  the  fourth  degree  polynomial  of 
Fediaevsky 's  method  are 

-f-=  I  e 

I  u, 


and 


»•(-  y  =  O 

-^=o  .+  y  =  rT 

Also,  from  the  energy  equation. 


-h 


3y 


=-2± 


(29) 


and  the  conditions,  v, 

7>y  ^  Iwy 
Moreover,  for^^ j  -  0 


-  0  and 


l-l-  0. 


it  follows  that  both 
O  akiocl  j  =  O  at-  y=0  . 


the  condition  at  y  -  6^  Is 

These  five  conditions  on  are  exactly  tha  same  as  on  -Pw  • 

Consequently,  if  and  .i  are  approximated  by  fourth  degree 

q  iw  "*>/  "X 

is  exactly  the  same  function  of  '>r  ns  is  of 

iw 


polynomials 


n.  Therefore,  the  relation  between  and  ’>p  is  assumed  to  be 

the  same  as  the  experimental  one  between  and  p.  Actually, 
however,  i  depends  slightly  on  the  Prandtl  number  (see  fig. 
7-55,  ref.'*(l)). 
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The  Prandtl  number  also  affects  the  ratio  in  equation 

(21).  Thus,  for  laminar  flow  over  a  flat  plate® 

-iL  rr  5Z3  ref.z) 

5* 

For  Pr  =*  .738,  the  value  used  in  the  present  analysis,  this 
relation  gives  1.106  for  .  For  turbulent  flow  ^  is 

S'  y 

probably  smaller  and  so,  as  an  allowable  approximation,  is 
taken  as  unity.  Equation  (21)  then  becomes 


Sj'R.j'Pr 


For  reasons  given  in  the  section,  "Friction  Coefficient  and 
Velocity  Profiles,"  (30)  is  written  as  , 

I  f  ~T 

5;T?es.?r  ^ 

The  values  y+^2  *^2  permitted  to  differ  from  y+,i  and  . 


To  find  the  ratio  of  the  eddy  thermal  diffusivity  to  the 
eddy  kinematic  viscosity,  4^  ,  the  pipe  flow  data  of  reference 

C»H 

(3)  were  used.  First,  the  values  of  given  in  reference  (3) 

for  five  Reynolds  numbers  at  each  of  twelve  values  of  y/r  were 
averaged.  This  was  allowable  because  there  was  no  consistent 
effect  of  Reynolds  number  on  at  the  twelve  stations.  A 

smooth  curve  of  against  y/r  was  then  drawn  through  these 

average  values  located  at  the  twelve  values  of  y/r.  Then  a 
curve  of  for  pipe  flow  was  obtained  by  first  averaging 

Laufer's  and  Runner's  data  given  in  figure  7-39  of  reference  (1) 
at  each  of  1  values  of  y/r  and  then  drawing  a  smooth  curve 
through  the  average  values.  Laufer's  and  Runner's  data  are  the 
same  from  the  wall  to  the  position  of  the  maximum  in  but 

Laufer's  data  are  slightly  higher  beyond  this  point.  Although 
Ref  is  estimated  to  be  about  6.3  x  10^  for  Laufer's  data  and 
about  6.6  X  10^  for  Runner's  data,  a  ratio  of  about  9.5,  it  is 
not  certain  that  this  is  the  reason  for  the  difference.  At 
each  value  of  y/6  or  y/r,  the  ratio  (4^ up)  then 

calculated  with  the  help  of  figure  2.  This  ratio  was  calculated 
because  the  eddy  kinematic  viscosity  approaches  zero  at  the 
outer  edge  of  a  boundary  lay^,  whereas,  it  does  not  at  the 
center  of  a  pipe.  The  ratio^^^^ was  plotted  against 

y/6  and  a  smooth  curve  drawn  througn  the  plotted  values.  The 
assumption  was  then  made  that  a  change  in  flow  from  pipe  to 
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plate  has  the  same  effect  on  e^j  as  it  does  on  Therefore,  r-^ 

was  obtained  by  multiplying;  the  ordinates  of  the  curve  for  ^  • 

obtained  from  the  pipe  flow  data  of  reference  (3)  by /£»  \ 

€k 

The  variation  of  with  y/6  for  plate  flow  was  then  calculated 

and  is  shown  in  figure  12.  For  n  S  •06,4^  Is  taken  as  unity. 

This  means  that  — is  assumed  to  be  £2(^1)  for  r|  <  ,06  with  the 
result  that  ^  , 

-!n?=-595n  CTi.o«) 

Consequently,  as  a  function  of  y^  has  a  horizontal  portion 


just  as 


Near  the  wall  the  proper  variable  is  y+  instead  of  ri.  For 

this  region  the  data  of  references  (6)  and  (7)  for  flow  of  air 
in  a  channel  were  used  to  obtain  the  connection  between  and 

y+.  The  values  of  ij^y  and  y+  were  calculated  by  using  p  ■■  22 

X  10“^  slugs  per  cubic  foot  and  v  -  1.8  x  10“^  square  feet  per 
second.  The  results  are  shown  in  figure  13  for  the  entire 
range  of  y.  from  the  wall  to  the  center  line  of  the  channel. 
Whore  Is  a  function  of  y^  the  data  fall  almost  on  fi,  except 

perhaps  between  about  3  and  17  where  the  ^  data  lie  somewhat 
below  fi .  The  agreement  between  these  ^  data  and  fi  is 

‘*•1^  6k 

believed  to  be  sufficiently  close  to  allow  fj^  to  represent 

as  well  as  .  This  means  that  ^  ••  1  near  the  wall  ” 

u»y  fern 

where  the  distance  variable  is  y^.  Instead  of  n. 

The  value  of  y^  2  corresponding  value  of  ri2  in  (31) 

ao'e  calculated  in  the  same  way  as  y^.  and  111.  To  do  this, 
figures  2  and  12  are  used  to  calculate  /ha.  )  \  as  a 

vlXvy  /  v  6  ipn  / 

function  of  y^  by  use  of  (9)  and  (11).  The  value  of  y^  at  which 
f  £hn  \  ^  equals  fi  is  then  found.  Below  this  value  of 

v»u»y/  '  <»«K/ 

y+  called  y+.2.  uTv  Is  represented  by  fi;  above  y+  2  figures 

2  and  12  are  used  to  got  (jjfyX  as  a  function  of  ti.  In 

Table  V  are  given  the  values  of  y^.  g  and  t\2  for  the  three  values 
of  Ref  for  which  y+^2  less  than’ 52. 


The  data  in  figure  13  indicate  that  the  assumption  that 
there  is  a  quantity  y+^1  or  y^.  2  seems  to  bo  correct.  The  method 
for  the  calculation  of’y+  i  or’y+  2  also  appears  correct.  That 
is,  the  data  show  that  below  the  Intersection  of  4^(1)  with 
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fl(y+),  depends  on  y^.  Above  the  intersection,  depends 

on  y^  and  on  Ref.  That  is,  above  the  intersection  the  variable 
is  and  it  depends  on  n.  As  expected,  there  is  a  smooth 

transition  between  if  and  f2  rather  than  the  corner  that 

results  in  the  present  method  for  finding  y+,1  and  y^  2- 
Ref  150  and  264,  the  experimental  data  leave  fi  cloie  to  the 
calculated  value  of  y+  2*  ^or  Ref  -  701,  there  is  an  appreciable 
difference.  Above  y^  2  calculated  curves  for  all  three 

values  of  Ref  are  below  the  experimental  data.  Because  the 
calculated  curves  are  for  flow  over  a  plate  whereas  the  points 
are  for  a  channel,  some  difference  is  to  be  expected.  Thus, 
although  the  value  of  -Siu  for  a  plate  approaches  zero  at  the 

outer  edge  of  the  boundary  layer,  the  va  jc  *  for  a  chan¬ 

nel  is  not  zero  at  the  center  line.  This  i.  y  the  dashed 
curves  in  figure  13  go  through  -  0  at  y^.  -  -  f,  but  the 

points  do  not.  For  y+  i  <  y+  <  Ref,  but  not  close  to  y+  -  Ref, 
only  a  small  part  of  ttie  difference  between  the  curves  and  the 
data  points  can  be  removed  by  multiplying  for  plate  flow  by 

the  reciprocal  of  the  ratio  \,Uvr// r  y  that  was  used  to 
convert  the  pipe  flow  ^data  of  reference  (3)  to  plate  flow 

data.  The  reason  for  the  remaining  difference  is  not  clear 

at  present. 

By  use  of  equation  (31),  the  function  fj,  figures  2  and  12, 
and  the  method  for  finding  y+,2»  the  local  Stanton  number,  St, 
was  calculated  for  the  five  values  of  Ref  listed  in  Table  I. 

The  results  are  given  in  Table  VI  and  also  shown  in  f Igure  9  as 
a  function  of  Req.  Also  given  in  Table  VI  and  figure  14  are  the 
values  of  the  ratio 

ZSr 

The  ratio  was  calculated  by  noting  that  from  the 

definition  of  Sx  and  of  Sf  it  follows  that 


5-r  = 


-  -T^cx 
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Upon  the  use  of  (19)  and  (21),  (33)  becomes 


SlIi  ■= 


For  0  <  Re(^  <  100,  the  ratio  was  expressed  in  an  analytic 

Cl 

form  by  assuming  that 

r  _  K I 

-n  •'I  (35) 

The  constants  kj  and  ni  were  determined  from  the  values  of  S-^ 
for  Reg  equal  to  27.7  and  191  (see  Table  VI)  and  were  found  to 
equal  .109  and  .601,  respectively.  When  the  results  for  kg  and 
nQ  are  used,  the  expression  for  ZSt.  for  Re^  <  100  becomes 

Z  S-tr  _  1*  "Sfi  ' _ 


,OZ2.S' 


9  V 

Equation  (36)  was  used  with  (34)  to  obtain  1.274  for 

Re0  100.  For  values  of  Rep  larger  than  100,  was  calcu¬ 

lated  by  numerical  integration  by  use  of  (34)  and  figure  (14). 
The  result  is  shown  in  figure  (14). 

Ft 

The  temperature  profiles  were  calculated  by  putting  "y  -  1, 
solving  (26)  for  ,  and  integrating  between  0  and  r|.  Thus, 

^  A  4, 


■t=^ 


w. 


_ 

1  +  45;  W' 


The  temperature  profiles  calculated  by  use  of  (38)  are  shown  in 
figure  15  and  tabulated  in  Table  VII.  Also  shown  is  the  tempera¬ 
ture  profile  for  laminar  flow.  This  profile  is  obtained  by 


letting 


i  in  (37). 


The  result  is 


or,  after  using  (17) 


•t  =  .S't'Reg-'Pr  I 


-L  = 
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From  (30)  with  from  (16)  it  follows 


^TtPr _ 


=  1 


^  (40) 

is  assumeri  to  be  nepligible. 


Therefore,  from  (39)  t  -  4*  .  That  is,  the  present  method 
predicts  that  for  laminar  flow  the  velocity  and  temperature 
profiles  are  the  same. 

Mass- Transfer  Coefficient  and  Concentration  Profile 

The  mass- transfer  coefficient  and  concentration  profile 
are  calculated  in  the  same  way  as  the  heat- transfer  coefficient 
and  temperature  profile.  For  a  binary  mixture  in  laminar  flow 
the  mass  transfer  by  diffusion  in  the  direction  normal  to  the 
wall  is,  _ 

W)  = 

where  the  contribution  of  t  "5^ 

For  turbulent  flow  it  is  assumed  that  an  eddy  dlffuslvlty  exists 
such  that 

=‘’j®  (41) 

After  introduction  of  the  Schmidt  number  ,  and  the  quantity 

~Z!  _c, 

equation  (41)  becomes  -«  /- 

In  the  same  way  as  in  the  section  concerning  heat  transfer,  it 
is  assumed  here  that  the  concentration  boundary  layer  thickness 
differs  from  the  velocity  boundary  layer  thickness  by  less  than 
10  percent  and  that  it  is  thus  permissible  to  ted^e  the  two 
thicknesses  equal. 

When  (42)  is  divided  through  by  mi  and  the  mass- transfer 
coefficient  ’ 


'(j. 


\y\ 


;omes 


C, 


—  m 


I  tj 


is  introduced,  the  result  is 

ywi  _  _ I 

■>^1. 


'w 


5b: 


2L 

n 


(43) 


(44) 


From  (44)  it  follows  that 

/*  I 


C-T?rs5'c  j  I 


(45) 
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Information  concerning:  i-1^  as  a  function  of  -n  is  obtained 
by  noting  that 


VWt  _  I 


n4-  y  =  o 


=  O  /=  S' 

vvt|  ^  » 

Moreover,  from  the  diffusion  equation 

[  r  ay 

it  follows  that  , 


Xi.UWQ  Uliau  k 

(^jw=0 

when  the  diffusion  velocity  is  small  enough  for  v^  ■■  0. 
If,  in  addition,  ©1-  0  and  ©i-  0,  then  also 


and  V  &  /  *■ 

Because  these  boundary  conditions  are  the  same  as  those  for  -r  , 

it  is  assumed  that  the  non-dimensional  concentration  distribution 
is  the  same  as  the  non-dimensional  heat- transfer  distribution. 

It  is  also  assumed  that  JX  -  1  and  that  Sg  -  Pr,  From  these 
assumptions  it  follows  that  _ 

Moreover,  the  non-dimensional  concentration  profile  is  then 
also  the  same  as  the  non-dimensional  temperature  profile.  If 
either  4^  /  1  or  S-^Pr,  the  method  used  to  calculate  the  heat 

Ni*  » 

transfer  coefficient  and  temperature  profile  can  be  used  to 
calculate  the  mass-transfer  coefficient  and  concentration  profile. 


DISCUSSION 

Friction  Coefficient  and  Velocity  Profile 

One  of  the  methods  for  obtaining  the  logarithmic  friction 
formula  and  velocity  profile  is  Millikan's  method  of  overlap. 
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In  this  method  It  is  shown  that  if  there  is  a  region  of  the 
boundary  layer  in  which 


^ 


are  both  valid,  tliej*  thf  iction  formula  is  the  logarithmic  one 
and  the  velocity  profile  is  logarithmic  in  the  region  of  overlap. 


In  the  present  analysis  the  logarithmic  friction  formula 
and  velocity  profile  also  follow  from  an  overlap  condition. 
Here,  however,  the  condition  is  on  the  eddy  kinematic  viscosity 
instead  of  on  velocity  profile  functions.  Thus,  to  obtain  a 
friction  formula  for  large  Reynolds  numbers  equation  (10)  is 
used.  For  large  Reynolds  numbers 


-rkf  ^ 

The  value  of  ni  is  thus  very  small.  Consequently,  Is  almost 

unity  lor  s  52.  Moreover,  the  indication  from  figure  1  is 
that  is  practically  unity  for  ri  less  than  about  .01,  In 
addition,  for  large  enough  values  of  Rey,  equation  (12)  applies 
for  a  range  of  n  that  begins  at  and  extends  outward  at 

least  as  far  as  .01.  Equation  (10)  then  becomes  | 

r ^ I  f  cjyf  f  ^  <jv; 

For  large  Re^  the  quantity  (t  -f  Reff2)  in  the  last  term  of  (48) 
is  almost  equal  to  Re^fo  except  very  close  to  -  1.  For 
example,  at  Ref  -  5  x  iQS,  Reff2  "  1965  at  ri  -  .01.  It  then 
increases  to  about  34,000  at  ri  -  .31  and  decreases  to  sero  at 
Ti  -  1;  at  Ti  -  .98,  however,  it  is  equal  to  850.  For  large 
Reynolds  numbers  the  last  term  in  (48)  can  be  written  as 

l-tv  ,  I 

y'— S  j 

.o| 

where  (1-^)  is  the  upper  limit  of  the  range  of  e  in  which  Reff2 
is  much  larger  than  unity.  As  Ref  Increases,  A  decreases  and 
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becomes  negligible.  For  large  Ref  the  last  term  in  (48)  is  then 
approximately  Independent  of  Re#.  Its  value  is  14.27.  The  first 
term  in  (48)  has  the  value  13.94.  When  these  two  numerical 
values  are  used  and  the  second  term  of  (48)  is  Integrated 
analytically,  the  result  is 


^  —  Z6.ZI  +  ^ 

For  large  Reynolds  numbers 

Ih  -1--0  I  K 'Rt^^  ~  In  (f-olKKe^:) 
Equation  (49)  then  becomes 


—  4-4>.31 


(K-=.-595j 


(49) 


(50) 


Note  that  the  logarithmic  form  requires  only  that  there  be  a 
region  of  the  boundary  layer,  no  matter  how  small,  In  which  (12) 
and  ^  •  1  are  both  valid. 

Equation  (50)  has  the  form  of  von  Karman's  friction  formula 
(ref.  (8)).  Although  (50)  was  derived  for  very  large  Reynolds 
numbers  the  value  of  Cf/2  obtained  from  It  differs  from  the  value 
given  in  Table  II  by  less  than  2  percent  for  Ref  2960. 

By  writing  (50)  in  the  form 

■=■  2.S^  io  iT? In  -4- 1 

y  _ 

and  using  the  value  of  'S’  tor  Ref  -  5  x  10®  from  Table  IV, 
there  is  obtained  the  relation 

1=^  =  2.54  U  r? (5X, 

The  value  of  Cf/2  obtained  from  (51)  differs  from  that  in  Table  II 
by  almost  7  percent  for  Ref  —  2960.  Equation  (51)  has  larger 
errors  than  (50)  at  small  Reynolds  numbers  because  the  ratio 
depends  on  Reynolds  number  (see  Table  IV)  and  the  value  of 
for  Ref  5  X  10®  was  used  to  get  (51).  Because  formulas  of  the 
logarithmic  type  are  Inconvenient  for  computation,  it  is  usually 
better  either  to  use  the  curves  of  figures  9  or  10  or  convenient 
approximations . 
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The  friction  coefficient  was  calculated  for  turbulent  flow 
over  the  entire  range  of  Reynolds  number  by  use  of  the  one  non- 
dimensional  shear  distribution  shown  in  figure  1,  The  good  agree¬ 
ment  between  the  friction  coefficient  calculated  by  the  present 
method  and  the  accepted  friction  coefficient  (see  figs.  5  and  10) 
indicates  that  the  non-dimensional  shear  distribution  is  approxi¬ 
mately  independent  of  Reynolds  number.  Moreover,  when  the  value 
of  ^  is  calculated  for  the  laminar  velocity  profile,  and 
substituted  into  (16),  the  result  is 

Z. 

This  expression  differs  from  the  exact  one 


W  _  .ZZOS 


by  less  than  2  1/2  percent.  Because  (52)  was  obtained  by  use  of 
the  same  shear  distribution  as  was  used  for  turbulent  flow,  the 
inference  is  that  the  non-dimensional  shear  distribution  on  a 


flat  plate  is  approximately  the  same  for  laminar  as  for  turbulent 
flow,  Fediaevsky's  method  (ref.  (4))  predicts  them  to  be  identical. 


In  contrast  to  the  usual  method  In  which  the  velocity  pro¬ 
files  are  calculated  by  patching  the  results  of  a  number  of 
approaches,  each  yielding  a  portion  of  the  velocity  profile,  the 
present  analysis  calculates  the  entire  velocity  profile  by  one 
approach.  Moreover,  the  present  method  yields  many  of  the 
results  that  are  presently  accepted.  For  example,  in  the  region 
very  close  to  the  wall  equation  (14)  applies.  Here  f.  «  1  and 
^  *  1.  Then  (14)  becomes 

^  I - -  r 

4'=  fF 

or 

“  j-h  (54) 


This  is  the  well-known  velocity  profile  of  the  laminar  sublayer 
(ref.  (1)). 

From  equations  (1)  and  (2)  it  follows  that  the  ratio  of  the 
turbulent  shear,  ,  to  the  laminar  shear,  ,  is  ^  , 

or  y+fj,.  The  value  of  unity  for  this  ratio  has  often  been  sug¬ 
gested  as  a  criterion  for  the  edge  of  the  laminar  sublayer.  By 
use  of  figure  3  it  is  found  that  this  ratio  is  equal  to  unity 
when  y+  is  about  9.9,  a  value  in  the  range  of  accepted  values 
for  the  outer  edge  of  the  laminar  sublayer. 

T 

For  somewhat  larger  values  of  y^.  the  ratio  tv,  is  still 
approximately  equal  to  unity  if  Re^  is  sufficiently  large. 
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Equation  (55)  is  the  law  of  the  wall  (ref.  (1)),  Its  upper 
limit  is  either  the  upper  limit  for  fj  or  for  T  'S' 1 ,  whichever 
value  of  y+  is  smaller. 


When  Re-  is  large  enough  for  a  region  of  overlap  to  exist, 
equation  (12)  is  valid  for  a  range  of  r|  whose  upper  limit  is 
near  .06  (see  figure  2).  Equation  (15)  then  becomes 


The  smallest  value  of  (1  +  K^Ref)  is  (1  +  52K) .  But  In  (1  +  52K) 
differs  from  In  52K  by  less  than  one  percent  for  K  -  .393. 
Consequently,  (57)  can  be  written  as 


—  15.  94  U  In 

f? 

or  as 


~h  Z.5-4 


-f-3.a  9 


Equation  (58)  is  the  logarithmic  profile. 


CA-i--"*-) 


22 


NOLTR  63-77 


Equation  (58)  was  obtained  by  requiring  that  *  1  and 

that  (12)  both  apply.  Consequently,  it  would  seem  that  the 
profile  should  not  be  logarithmic  for  ti  ,06  or  so;  n  -  .06  is 
roughly  the  upper  limit  for  (12).  When,  however,  (9)  is  used 
with  figure  7  it  is  found  that  the  velocity  profiles  do  not 
differ  noticeably  from  the  logarithmic  profile  unless  ri  is  much 
larger  than  .06.  For  example,  the  values  of  t|  for  departure  from 
the  logarithmic  profile  vary  from  about  ,17  for  Ref  -  2960  to 
about  .20  for  Ref  -  5  x  10^,  For  Ref  equal  to  200  there  appears 
to  be  very  little,  if  any,  logarithmic  portion  in  figure  7;  the 
lower  values  of  Rof  definitely  have  none.  The  velocity  profiles 
thus  have  a  logarithmic  appearing  portion  for  values  of  t)  larger 
than  a  strict  application  of  the  requirements  indicate.  The 
reason  for  this  behavior  is  probably  the  slow  decrease  of  Sen- 
from  Kn  (see  fig.  2)  combined  with  the  slow  decrease  of  X 
from  unity  for  n  less  than  about  .2  (see  fig.  1).  ^ 

As  the  Reynolds  number  decreases,  the  value  of  ni,  which 
for  large  Reynolds  numbers  is  equal  to  ^  ,  Increases.  The 

increase  in  rii  decreases  the  range  of  ti  ?*(^•r  which  (12)  is  valid 
and,  as  a  result,  decreases  the  extent  of  the  logarithmic  portion 
of  the  velocity  profile.  If  the  upper  limit  of  (12)  is  taken  at 
Ti  ~  .06  (see  fig.  2),  there  will  be  no  logarithmic  portion  of  the 
velocity  profile  when  Re^  becomes  less  than  867.  At  this  value 
of  Ref,  rii  exceeds  ,06.  The  value  867  for  Ref  corresponds  to  a 
value  of  2000  for  Req.  Note  that  the  lower  limit  for  the  loga¬ 
rithmic  velocity  profile  is  also  the  lower  limit  for  the  loga¬ 
rithmic  friction  formula,  equation  (50) , 

Although,  strictly  speaking,  there  is  no  logarithmic  portion 
of  the  velocity  profile  for  Ref  <  867,  the  discussion  concerning 
the  extent  of  the  logarithmic  profile  indicated  that  according 
to  flg;ure  7  a  velocity  profile  has  a  logarithmic  appearing 
portion  down  to  Ref  In  the  neighborhood  of  200.  For  Ref  -  200, 
Req  is  about  450.  Preston  (ref.  (9))  took  the  lower  edge  of  the 
logarithmic  portion  of  the  profile  at  y+  -  30  and  the  upper  edge 
at  n  -  .20  and  found  that  the  extent  of  the  logarithmic  portion 
shrank  to  zero  at  Rcf  -  150.  From  this  value,  Preston  obtained 
389  for  the  value  of  Req  below  which  there  is  no  logarithmic 
portion  of  the  velocity  profile.  From  figure  3  It  seems  that  30 
is  too  small  a  value  of  y^  for  the  logarithmic  profile  to  hold. 
Here  again,  however,  the  velocity  profile  departs  slowly  from 
the  logarithmic  type  (see  also  ref.  (10))  so  that  it  has  the 
logarithmic  appearance  for  smaller  values  of  y+  than  is  to  be 
expected . 

Beyond  the  logarithmic  portion  of  the  velocity  profile  lies 
the  region  of  the  velocity  defect  formula  (ref.  (1)).  To  obtain 
this  formula  from  the  present  analysis,  equation  (15)  is  written 
for  n  ”  1 ,  thus. 
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When  (15)  is  subtracted  from  (59),  the  result  is 


l-4> 

W 


-  '~^Cr 


Although  and  <2  taken  to  be  functions  of  n  alone,  the 

indication  from  (60)  is  that-p»  depends  on  Rsf  as  well  as  on  ^ 

unless  Ref  is  large.  This  result  is  illustrated  by  the  profiles 
for  Rsf  "  25,  100,  and  200  in  figure  8, 

For  large  Rsf  the  quantity  (1  -f  Reff2)  is  approximately 
equal  to  Reff2  for  all  >  ’ll  except  at  t|  -  1,  where  f2  "  0. 
Therefore,  for  large  Ref  (60)  becomes  j 

As  indicated  by  the  discussion  following  equation  (48) ,  the  last 
term  of  (61)  becomes  negligible  for  large  Re^.  Consequently, 
for  large  Ref  (61)  becomes 

“  -A  ^*1  ) 

or,  in  the  usual  form. 


Equation  (61)  is  the  velocity  defect  formula  (ref.  (1)).  Its 
derivation  indicates  it  to  be  valid  only  for  sufficiently  large 
Reynolds  numbers.  The  indication  from  figure  8  is  that  (62)  is 
valid  for  the  two  largest  values  of  Re^. 
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The  present  analysis  yields  the  velocity  defect  formula  and 
the  law  of  the  wall.  Explicit  expressions  for  two  portions  of 
the  law  of  the  wall,  namely,  the  laminar  sublayer  and  the  loga¬ 
rithmic  portion  are  obtained.  Both  the  velocity  defect  and  the 
logarithmic  formula  require  a  sufficiently  large  Reynolds  number 
and,  in  addition,  the  logarithmic  formula  requires  that  be 

almost  unity. 


It  is  remarked  that  although  is  the  same  function  of  y^. 


for  both  pipe  and  plate  flow  and 


u.r 


for  rt  less  than  about  .2,  the  shear  ratio  is  not. 
distribution  in  a  pipe  is 

r.  ^  i-JL 

Tw  '  r 


is  the  same  function  of  rt 
The  shear 


rather  than  the  distribution  shown  In  figure  1.  Consequently, 
the  law  of  the  wall,  which  was  obtained  by  requiring  that  ^ 

« Wk/ 

be  almost  unity,  should  not  extend  out  as  far  from  the  wall  for 
pine  flow  as  for  plate  flow  (see  p,  617  of  ref.  (D).  Moreover, 
for  most  of  the  range  of  n  for  which  the  velocity  defect  law 
holds,  both  is;  and  S.  are  different  for  pipe  thp.n  for  plate  flow. 

Thus,  in  a  pipe  remains  equal  almost  to  its  maximum  value 

out  to  large  values  of  n  instead  of  dropping  to  zero. 
Consequently,  the  velocity  defect  formula  should  be  different 
for  pipe  than  for  plate  flow  (ref.  (1)). 


Heat- Transfer  Coefficient  and  Temperature  Profile 

The  ratios 
Rea,  the  ratio 


iSt 


3  ^ —  and  cp  are  shown  in  figure  14. 
1*1^  is  larger  than  because  i 


For  equal 
is  an  average 


of  (see  eq.  (34))  and  thus  includes  the  larger  values  of 


that  are  present  at  low  values  of  Reo. 
is  the  line 

2^  — 

Cr  Cp 


Also  shown  in  figure  14 

_a. 

T 


the  commonly  accepted  ratio  (p,  497,  ref.  (2)).  The  present 
analysis  results  in  a  ratio  or  that  is  larger  than 

Pr"2/3  low  Reynolds  numbers  and  smaller  at  high  Reynolds 

numbers.  The  behavior  at  very  large  Reynolds  numbers  is  obtained 
from  (31).  For  large  Reynolds  numbers  (31)  becomes 


1 


■S-A 


( 


~h 


(63) 
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q  Sz. 

or,  after  taking-^  -  1  for  <  .01  and  ±2  -  Kti  for  <  n  <  .01, 

equation  (63)  becomes  ,  ,  ” 


(64) 


J  I  f  y/,  ?r  I  4-K^-1?,^Fr  J  1+ 


For  Pr  -  .738,  calculations  gave 
The  last  Integral  Is  equal  to 


JS.6A 

A\  I  I 


for  the  first  integral. 


.  J, 

445.  ^ 


which,  for  large  Re^  is  equal  to  ror  t-r  -  ./aa.  me 

second  integral  is  Integrated  analytically.  Equation  (64)  then 
becomes 


1-^  n-if 


for  Pr  -  .738.  The 


S^'T^es"  L.  ^ 

or,  for  .OIK  Re  Pr  »  1, 


(^1  +.  ol  KFp^Fr}— 


I  r  -| 

.^.b4  lv^'~'5er  (65) 

When  (50)  is  used  for j|Cj/2,'  the  result  is 

■ZSt  i’.s4  ^  *  _ 

— - —  - ^  (66) 

-h  3.0^ 

For  very  large  Reynolds  numbers  (66)  becomes 

=  (  (67) 

The  result  that  the  ratio  approaches  unity  does  not  depend 

on  the  value  of  K  or  of  Pr.  Because  of  the  logarithm  in  the 

numerator  and  denominator  of  (66) ,  the  ratio  approaches  unity 

very  slowly.  For  example,  for  Re«  -  lOlO  zit  -  1.053. 

* 

From  equations  (13)  and  (37)  it  follows  that  if  Pr  -  1 

for  all  n.  then  ^ 

4-  J2  5*.?;- 

Ip 

or,  because  t—^)—  latn—  1 

-  1 
"  ‘ 


for  all  r|«  then 
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This  also  follows  from  (5)  and  (30).  Therefore,  the  velocity 
and  temperature  profiles  are  Identical  when  Pr  -  1  for  all  ri. 
Because  depends  on  this  cannot  be  so.  Cc%sequently ,  even  if 
-  ^""and  Pr  -  1  the  velocity  and  temperature  profiles  are 
not  Identical. 

i  ^1^  n  j 

t  ^  “  O,  and 


(68) 

The  correct  result  (ref.  (2))  is  that 

(69) 


For  laminar  flow,  it  follows  from 
from  (16)  and  (30)  that 

2St  _  _L_ 
~cr7  — Pr- 


A  comparison  of  equation  (16)  with  equation  (30)  for  "tT  -  0, 
indicates  that  the  cause  of  the  error  is  the  assumption  that 
'  "fli  should  depend  on  Pr.  Thus,  it  can  be  shown 

from  pages  120  and  313  of  reference  (2)  that  for  laminar  flat- 
plate  flow  _ 


( 

,■532. 


o 


and  that 


(70) 


The  quantity  Pr“^/^  in  (70)  accounts  for  the  difference  between 
(68)  and  (69).  The  inference  is  therefore  that  the  effect  of 
Pr  on  the  ratio  in  the  present  analysis  is  Inexact.  When 

the  Prandtl  number^  is  near  unity  the  error  is  small.  Thus,  for 
Pr  -  .738,  the  error  in  (68)  is  about  12  percent. 


Results  similar  to  those  obtained  for  the  velocity  profiles 
follow  for  the  temperature  profiles.  Thus,  very  close  to  the 
wall,  fi  «  1  and  ^  ^  1.  Equation  (38)  then  becomes 
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Define  the  quantity  t«  as 


■t* - 


(see  ref.  (11), p.  823) 


Equation  (71)  then  becomes 


where 

This  Is  the  linear  portion  of  the  temperature  profile  very  close 
to  the  wall. 


-t  ,  “ 


The  turbulent  shear  is  equal  to  the  laminar  shear  where 
y+fj^  Is  equal  to  unity.  This  occurs  at  about  y.  -  9.9.  From 
the  denominator  of  the  first  Integral  of  (38)  It  follows  that 
the  turbulent  heat-transfer  Is  equal  to  the  laminar  where 
y^fiPr  equals  unity.  For  Pr  <  1,  the  so-called  edge  of  the 
laminar  temperature  layer  Is  larger  than  the  edge  of  the  laminar 
velocity  layer.  Thus,  for  Pr  -•  .738  the  edge  of  the  temperature 
laminar  sublayer  Is  at  about  V-t- **  11.4,  a  valvte  slightly 
larger  than  the  value  of  about  9.9  for  the  velocity  layer. 

For  y^  too  large  for  (72)  but  small  enough  for  to  be 
approximately  equal  to  unity,  (38)  becomes 

■t  = 

/IF  Jo 


then 


Equation  (73)  Is  the  law  of  the  wall  for  the  temperature  pro¬ 
file.  Like  the  law  of  the  wall  for  the  velocity  profile.  Its 
upper  limit  Is  either  the  upper  limit  for  fj  or  f or  ^  -  1, 
whichever  is  smaller. 

For  the  range  of  rt  for  which  (12)  is  valid,  (38)  can  be 
written  as 


StTPr 


-U-R- 


"Ro  P  I  - - - 7 — 
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For  large  enough  Ref,  Is  almost  unity  up  to  values  of  n 

large  enough  to  include  the  upper  limit  for  (12) ,  which  is  at 
about  .06.  Figure  12  Indicates  that  for  -n  between  zero  and 
about  .06,  ^  is  almost  unity.  Therefore,  (74)  can  be  written 
as  ^ 


ex 


“b  = 


=:  d^H 


■h  S^l^e/Fr 


The  first  integral  in  (75)  has  the  value  15.88  for  Pr  -  .738. 
Equation  (75)  then  becomes  ^ 


-t=  unt-^ 

f¥  k/F 


I  4-  K 


or 


-Tiv 


=  \UI  +--k-  L  [l  U[|+S^K^.^3»3 


or,  with  K  ••  .393  and  Pr  ■  .738, 


"t*  ,  =  7.  S'4lvnri  4-  Z  54  Iv^'^er  -f- (.So  /ygx 

-  *1  i .  oC) 

Equation  (76)  is  the  logarithmic  temperature  profile.  ^ 


A  temperature  defect  formula  can  be  obtained  by  writing 
(38)  for  n  -  1,  thus 

I  —  I 


i. 


When  (38)  Is  subtracted  from  (77) ,  the  result  is 


I— 


‘1 


or 


-  w  Tl.  - ^ 
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or,  after  using  the  definitions  of  t,  t^,  and  Ref, 

We-T 

The  indication  from  (78)  is  that*~^;jj — depends  on  Ref  and  Pr  as 
well  as  on  n.  For  large  Ref,  however,  (78)  becomes 

"Te-T  f  ~ -  , 


j  H  j  I  \  (79) 

where  (1-A)  is  the  upper  limit  of  the  range  of  n  in  which 

^  Pr  is  much  larger  than  unity.  For  very  large  Ref,  the 

last  term  In  (79)  becomes  negligible  and  (79)  becomes 

Equation  (80)  is  a  temperature  defect  profile.  The  velocity 
defect  profile  (eq.  (6^))  requires  only  a  sufficiently  large 
Ref.-  In  addition  to  this  requirement,  (80)  requires  that  J- 

and  both  depend  on  n  alone.  Because  Is  known  to  depend 
on  Pr  and  because  ^  probably  also  does  (see  p.  552,  ref.  (1)), 


on  Pr  and  because  ^  probably  also  does  (see  p.  55J 
equation  (80)  should’*'more  properly  be  written  as 

^=4  (VP.; 

Mass- Transfer  Coefficient  and  Concentration  Profile 


The  case  treated  corresponds  to  a  turbulent  boundary  layer 
composed  of  two  species  flowing  over  a  plate  on  which  one  of  the 
species,  say  Sf,  is  condensing  or  evaporating.  The  concentration 
of  Sf  is  supposed  to  be  so  small  that  the  velocity  at  the  wall 
caused  by  the  evaporation  or  condensation  of  Sj^  is  not  large 
enough  to  make  the  concentration  profile  differ  from  the  shear 
or  heat- transfer  profile  for  v^  -  0.  If  C^  is  the  concentration 
of  the  evaporating  or  condensing  species,  Vw  is  given  by 
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All  the  results  obtained  for  the  temperature  profiles  and 
heat- transfer  coefficients  can  be  converted  to  the  same  results 
for  the  concentration  profiles  and  raass-transfer  coefficients  by 
making  the  substitutions 

*t^6c 

t  — ? 

Pr-»S^ 

c 

«h-*®d 

q— ^  mi 

k  — «»  pD 

T— -Cl 


It  has  been  shown  that  the  friction  coefficient  for  a  con¬ 
stant  property,  zero  pressure  gradient,  turbulent  boundary  layer 
can  be  calculated  over  the  entire  range  of  Reynolds  numbers  by 
one  method.  The  method  is  to  use  the  relation  between  the  local 
shear,  the  local  sum  of  the  molecular  and  eddy  viscosity,  and 
the  local  velocity  gradient.  The  eddy  viscosity  across  the 
boundary  layer  is  obtained  by  Joining  a  distribution  for  the 
wall  region  to  one  for  the  outer  region.  The  same  procedure 
yields  the  velocity  profile  from  the  wall  to  the  outer  edge  of 
the  boundary  layer  without  using  the  concept  of  laminar  sublayer, 
transition  region,  logarithmic  region,  etc. 

The  heat  and  mass-transfer  coefficient  and  the  temperature 
and  concentration  profile  are  calculated  by  a  similar  method. 
More  approximations,  however,  are  used  than  to  calculate  the 
friction  coefficient  and  velocity  profile.  Therefore,  these 
results  are  probably  not  as  accurate  as  those  for  the  velocity 
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boundary  layer.  For  example,  although  the  calculated  friction 
coefficient,  Cf/2,  agrees  closely  with  accepted  values,  the 
ratio  differs  slightly  from  the  accepted  value,  ,  and 

decreases  with  increasing  Reynolds  number. 
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FIG.  I  RATIO  OF  LOCAL  TO  WALL  SHEAR  AS  A  FUNCTION  OF  RATIO  OF  WALL 
DISTANCE  TO  BOUNDARY  LAYER  THICKNESS 


FIG.  2  RECIPROCAL  OF  EDDY  KINEMATIC  VISCOSITY  REYNOLDS  NUMBER  AS  A  FUNCTION 
OF  RATIO  OF  WALL  DISTANCE  TO  BOUNDARY  LAYER  THICKNESS 
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FIG.  3  RECIPROCAL  OF  EDDY  KINEMATIC  VISCOSITY  WALL  DISTANCE  REYNOLDS  NUMBER  AS  A  FUNCTION 
OF  THE  KINEMATIC  VISCOSITY  WALL  DISTANCE  REYNOLDS  NUMBER 
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FIG. 5  LOCAL  FRICTION  COEFFICIENT  AS  A  FUNCTION  OF  BOUNDARY  LAYER 
FRICTION  REYNOLDS  NUMBER 
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REYNOLDS  NUMBERS 


▼ 


I 


FIG.7A  RATIO  OF  LOCAL  VELOCITY  TO  FRICTION  VELOCITY  AS  A  FUNCTION  OF  WALL  DISTANCE 
REYNOLDS  NUMBER  FOR  SIX  FRICTION  REYNOLDS  NUMBERS 
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FIG.7B 
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6.8  VELOCITY  DEFECT  AS  A  FUNCTION  OF  RATIO  OF  WALL  DISTANCE  TO  BOUNDARY- LAYER 
THICKNESS  FOR  SIX  FRICTION  REYNOLDS  NUMBERS 
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FIG.  9  LOCAL  FRICTION  COEFFICIENT  AND  LOCAL  STANTON  NUMBER  AS  A  FUNCTION  OF 
BOUNDARY  LAYER  MOMENTUM  THICKNESS  REYNOLDS  NUMBER 
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FIG.  10  AVERAGE  FRICTION  COEFFICIENT  AS  A  FUNCTION  OF  REYNOLDS  NUMBER  BASED  ON 
DISTANCE  TO  LEADING  EDGE 
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BOUNDARY  LAYER  MOMENTUM  THICKNESS  REYNOLDS  NUMBER  AS  A  FUNCTION  OF 
REYNOLDS  NUMBER  BASED  ON  DISTANCE  TO  LEADING  EDGE 
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THICKNESS  REYNOLDS  NUMBER. 
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FIG.  15  NON-DIMENSIONAL  TEMPERATURE  PROFILES  FOR  SIX 
FRICTION  REYNOLDS  NUMBERS 
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TABLE  I 


Limit  of  Wall  ReRion  and 

Associated 

Eddy  Kinematic 

Viscosity 

Reynolds  Number  For 

Three  Friction  Reynolds  Numbers 

Re^ 

/  Em  \ 

1 

u*6 

25 

.127 

12.2 

.0615 

.487 

100 

.253 

26.6 

.0672 

.266 

200 

.318 

37.6 

.0598 

.188 
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Friction 

TABLE  II 

Coefficient  and  Boundary  Layer  Reynolds 

Number 

for  Six  Friction  Reynolds  Numbers 

Cf/2 

Reg 

25 

.0116 

2.32  X  102 

100 

.00378 

1.62  X  103 

200 

.00280 

3.78  X  103 

2,960 

.00143 

7.81  X  104 

22,400 

.000987 

7.13  X  105 

500,000 

.000633 

1.99  X  107 

Calculated  Velocity  Profiles 


NOLTR  63-77 


NOLTR  63-77 


TABLE  IV 


Computed  Values  ,  and 

l^or  Six  Friction  keynolds  Mumpers 


Ref 

r 

T 

5* 

"5" 

Re0  X  10”^ 

25 

.119 

.269 

2.25 

.0277 

100 

.118 

.211 

1.79 

.191 

200 

.117 

.189 

1.61 

.444 

2,960 

.102 

.138 

1.35 

7.98 

22,400 

.0922 

.119 

1.29 

65.7 

500 , 000 

.0791 

.0954 

1.21 

1571. 
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TABLE  V 


Limit  of  Wall  ReRJ on  and  Associated  Eddy  Thermal 
Diftuslvlty  Reynolds  Kumber  for  Three  friction 
Reynolds  jTumbers 


RCf 

i^) 

y+.2 

^^2 

25 

.140 

13.2 

.0736 

.528 

100 

.266 

28.4 

.0750 

.284 

200 

.328 

39.0 

.0639 

.195 
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TABLE  VI 

Local  Stanton  Number,  Ratio  of  Stanton  Number  to  Friction 
coefficient,  and  Boundary-Layer  Momentum  Thickness 
Reynolds  itumber  for  ^Ix  Friction  Reynolds  Numbers 


Rejf 

St 

25 

.0148 

100 

.00463 

200 

.00339 

2,960 

.00164 

22,400 

.00110 

500,000 

.00069 

^‘^t  Ref,  X  10~3 


Cf 

1.23 

.0277 

1.22 

.191 

1.21 

.444 

1.14 

7.98 

1.12 

65.7 

1.09 

1571. 

Calculated  Temperature  Profiles 
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